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Abstract. 

The nonlinear absorption of laser radiation of relativistic intensities in the 
underdense plasma by a mechanism of stimulated brc;msstrahhmg of electrons on the 
ions/nuclei is investigated in the low frequency approximation. Coefficient of nonlinear 
inverse-bremsstrahlung absorption is studied for relativistic Maxwellian plasma at 
asymptotically large values of laser fields and high temperatures of electrons. 
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1. Introduction 

In the last two decades, laser technologies have made a giant leap forward such that laser 
sources of ultrarelativistic intensity level can be attained |l] . Near-infrared laser beams 
are available up to intensities of 10^^ W/cm and much stronger lasers will be available 
in near future [2]. For such lasers the dimensionless relativistic invariant parameter 
of intensity = eEo/mcu) >> 1 (e -elementary charge, m - electron mass, Eq, u - 
electric field amplitude and frequency of a laser radiation, c -light speed in vacuum). 
The latter represents the work of the field on the one wavelength in the units of the 
particle rest energy. Interaction of such lasers with the matter at extreme conditions 
in ultrashort space-time scales have attracted broad interest over the last few years 
conditioned by a number of important apphcations, such as generation and probing 
of highenergy-density plasma [s] , ions acceleration and inertial confinement fusion |4] , 
vacuum nonlinear optics js], compact laser-plasma accelerators |6], etc. Generally, the 
interaction of such fields with the electrons in the presence of a third body makes 
available the revelation of many nonlinear relativistic electrodynamic phenomena. As a 
third body can serve ion and in the superintense laser fields one can observe relativistic 
above threshold ionization [T] and high order harmonic generation js] , electron-positron 
pairs production on nuclei |9], and multiphoton stimulated hremsstrahlung (SB) of 
electrons on the ions/nuclei [lO]. The latter is one of the fundamental processes at the 
interaction of superstrong laser pulses with plasma and under the some circumstances 
inverse-bremsstrahlung absorption may become dominant mechanism of absorption of 
strong electromagnetic (EM) radiation in underdense plasma. 

With the advent of lasers many pioneering papers have been devoted to the 
theoretical investigation of the electron-ion scattering processes in gas or plasma in 



the presence of a laser field using nonrelativistic 11 20 as well as relativistic 21 25 



considerations. The appearance of superpower ultrashort laser pulses of relativistic 



intensities has initiated new interest in SB in relativistic domain 26 -28 , where 



investigations were carried out mainly in the Born approximation over the scattering 



potential. Meanwhile for ions with the large charge and for the clusters 29 , when 
electron interaction with the entire dense cluster ion core that composed of a large 
number of ions is dominant, the Born approximation is not applicable. The theoretical 
description of SB in superstrong EM fields and scattering centers of large charges requires 
one to go beyond the scope of Born approximation over the scattering potential and 
the perturbation theory over laser field. In this context, when quantum effects are 



considerable, one can apply eikonal 23 30 or generalized eikonal approximation 31 



For the infrared and optical lasers, in the multiphoton interaction regime, one can apply 
classical theory and the main approximation in the classical theory is low frequency 



(LF) or impact approximation 16,17. LF approximation have been generalized for 
relativistic case in Refs. 24 ,25 , where the effect of an intense EM wave on the dynamics 
of SB and non-linear absorption of intense laser radiation by a monochromatic electron 
beam due to SB have been carried out. Regarding the absorption of an EM of relativistic 



Nonlinear stimulated hremsstrahlung 



3 



intensities in plasma due to inverse SB there are not investigations which could help to 
clarify the behavior of nonlinear absorption of an intense EM radiation at very large ^o- 
Hence, it is of interest to clear up how nonlinear SB effect will proceeds in the plasma, 
taking also into account initial relativism of plasma electrons. 

In the present paper inverse-bremsstrahlung absorption of an intense laser radiation 
in relativistic Maxwellian plasma is considered in the relativistic LF approximation 



24 , 25 . In this approximation one can consider as superstrong laser fields as well as 
scattering centers with large charges. The radiation power absorption in such plasma 
is investigated as for circularly polarized wave (CPW), as well as for linearly polarized 
wave (LPW). We consider the dependence of the absorption coefficient on the intensity 
and polarization of the laser radiation, as well as on the temperature of the plasma 
electrons. 

The organization of the paper is as follows: In Sec. II the relativistic absorption 
coefficient of the EM wave of arbitrary polarization and intensity due to the mechanism 
of SB process is presented. In Sec. Ill we consider the problem numerically along with 
derivation of asymptotic formulas for absorption coefficient. Conclusions are given in 
Sec. IV. 

2. Nonlinear inverse-bremsstrahlung absorption coefficient 

The absorption coefficient a for an EM radiation field of arbitrary intensity and 
polarization, in general case of the homogeneous ensemble of electrons of concentration 
Tie, with the arbitrary distribution function /(p) over momenta p, at the inverse 
hremsstrahlung on the scattering centers with concentration n^, can be represented 
in the form: 

71 f 

a = ^ / rfpo/(Po)W^, (1) 



/ „ 

where W is the classical energy absorbed by a single electron per unit time from the EM 
wave of intensity / due to SB process on the scattering centers. For the homogeneous 
scattering centers W ^ rii. For the generality, we assume Maxwellian plasma with the 
relativistic distribution function: 

g(po) ' 



eXp , ^^rp^ 

where ks is the Boltzmann's constant, is the temperature of electrons in plasma, £{po) 
is relativistic energy-momentum dispersion law of electrons, K2{x) is the McDonald's 
function; /(po) is normalized as 

f{po)d^Po = 1. (3) 

To obtain W for the SB process, the electron interaction with the scattering 
potential and EM wave in the LF approximation can be considered as independently 



proceeding processes 17,25 , separated into the following three stages, schematically 



Nonlinear stimulated hremsstrahlung 



4 



depicted in Fig. 1. Field free electron with energy 8q and momentum po interacts with 
the EM wave. The exact solution of the relativistic equation of motion of an electron in 



a plane EM wave is well known 10 . The energy and momentum in a plane EM wave 
field can be written as: 

P±W = P±(V^o) -e , (4) 

i/p(V') = vMo) + ^ I \ ^ TT^ 

2c[£{i)Q) - cz/p(?/'o)) 

X [e^ (A(^o) - A(^))^ - 2ecp(^o) (A(^o) - A(^))] , (5) 
^(^) = ^(^o) + cz/(pW-p(^o)), (6) 



where 



A(V^) = Ao(^)(eicosV^ + e2CsinV^) (7) 

is the vector potential of the EM wave of currier frequency u and slowly varying 
amplitude Aq{iP). Here ip = lot is the phase, r = t — i/r/c , z/ is an unit vector in 
the EM wave propagation direction, ei 2 are the unit polarization vectors, and arctan^ 
is the polarization angle. At the second stage the elastic scattering of the electron in 
the potential field takes place at the arbitrary, but certain phase ips of the EM wave. 
Thus, taking into account adiabatic turn on of the wave (A(^o) = 0) from Eqs. Q-Q 
before the scattering one can write 

P±Ws) = Po± + — - — , (8) 

up{4js) = ^Po + ^ [e^^\4^s) + 2cepoA(^,)] , (9) 

£(^,)=^o + cz/(p(V^.)-po), (10) 

where 

A = £{i)s) - cup{ijjs) = £o - ci^Po (11) 

is the integral of motion for a charged particle in the field of a plane EM wave. The 
mean energy of an electron in the wave-field before the scattering will be 

(^(^)).=^o + 4a^- (12) 
Then it takes place elastic scattering on the scattering center U (r). Due to the 
instantaneous interaction of the electron with the scattering potential the wave phase 
does not change its value during the scattering. The electron with initial momentum 
p{ips) is acquired the momentum p'{ips) {£'{i^s) = £{i^s)) after the scattering directly at 
the same phase ips of the wave, which can be defined from the generalized consideration 




of the elastic scattering. Thus, measuring the scattering angle from a direction p [tpg 
with corresponding azimuthal angle (p for the scattered momentum one can write 



p'x ii^s) 




sin -d cos (p 


Py ii^s) 


= pMR 


sin sin ip 


_ P'z ii^s) _ 




cost? 



(13) 



where R = Rz (v^o) Ry {1^0), Ry {1^0) and Rz (v^o) are basic rotation matrices about the y 
and z axes: 



R 



cos -^0 COS ipo — sin (po sin {}q cos 990 
cos "Sq sin y^o cos y^o sin {}q sin 
— sin -i^o cos -do 



(14) 



As a axis we take wave propagation direction u, 6q is the polar angle and (po is the 
azimuthal angle in the wave-polarization plane. 

At the third stage the electron again interacts only with the wave, moving in the 
wave field with the momentum and energy defined from Eqs. (|4])-(|6]): 

A(^,) - A(^) 



P±fW = p' i^s) -e- 

upfW = up'i^,) + ^ [e' (A(^,) - A{^)y 
-2cep'(^,) (A(^,) - A(^))] , 
£/(^) = £(V',) + cz/(p(^)-p'(^,)), 



(15) 

(16) 
(17) 



where 



A' = £f{^) - cupfiilj) = SiiPs) - cup'iilj,). (18) 
The mean energy of an electron in the wave field after the scattering will be 

X [e' {A\^,) + (A^)) - 2cep'(V^,)A(^,)] . (19) 

The energy change due to SB can be calculated as a difference of mean energy in the 
field before and after the scattering: 

AS i^,^,i;,,po) = {£fm~{£m^. 
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Taking into account Eqs. (12) and (19) we obtain: 

,A\i,,) + {A^) [ I 1 



A' A 



ecp'(V's)A(V^ ^ ecp{ips)A{iljs 



For the energy absorbed by a single electron per unit time from the EM wave due to 
SB process on the scattering centers one can write 



W 



ni_ 
2tt 



2-K 



(21) 



where v (■i/'s) = c^p{il)f^) / £{'il)s) is the velocity of an electron in the wave-field, p{ips) = 
^yS'^{ilJs) — m'^c'^/c is the momentum, and da{'d,p{iljs)) is the differential cross section 
of the elastic scattering in the potential field U (r). Taking into account that the main 
contribution in the integral (21) comes from the small angle scatterings one can write 



w 



27r 



dtps / 



d^A£ 



datr{'d,p{ips)), 



(22) 



where 

da,,{^, p{^,)) = (1 - cos ^) da{^, p(V's)) (23) 

is the transport differential cross section. For the Coulomb scattering centers with 
potential energy 

f/(r) = ^ 

r 

of electron interaction with ion of charge Ze, one can use relativistic cross section for 
elastic scattering at small angles 32 and make integration over d and to obtain: 



dip, 



X {£iiPs)A - m^c*) + cepiiPs)A{iP,)] 



-L 



cb; 



where 



^cb 



In 



v^ {tl)s)p{ips 



(24) 



(25) 



is the Coulomb logarithm. The latter has been obtained taking pmin = Ze^ /vp as a 
lower limit of the impact parameter, while for the upper limit we assume pmax = v/w. 
Taking into account Eqs. ([T]), ([7]), and (24) for the absorption coefficient we obtain: 



a 



X 



dpof (Po 



2n 



dips 



, A'{^s) + {A') 
2m2c4 



(^(^,)A 



2 4 

m c 



) + cep{ips)A{ip, 



m^c^Eiips 
K^p^iips) 



(26) 
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where 

/ = (1 + C'^)uj'^ All Stic. 



Thus, Eq. ( 26 ) represents nonhnear inverse-bremsstrahlung absorption coefficient a for 
an EM radiation field of arbitrary intensity and polarization, for homogeneous ensemble 
of electrons of concentration ng, with the arbitrary distribution function /(po) over 
momenta po- 

Note that the LP approximation in the intense laser field is applicable when 

A > As, (27) 

where A is the laser radiation wavelength and Ai^ = ^A;{,Te/47rnee^Z is the Debye 
screening length: 



A„|c,„1^7.43xlO=x,/^I|XL. (28) 



Besides, for an underdense plasma one should take into account condition u > Up, where 
Up = A'mif.e^ /m^ is the plasma frequency with "effective mass" of the relativistic 
electron in the EM wave [331: 



m* = m^/lTWiM- (29) 

Thus for CPW and for large one can write 

rie < m^ou^/Ane^ = 1.1 x 10^^ x x A-^(/im). (30) 

There is also limitation on the pulse duration r of an EM wave. The SB should be the 
main mechanism, which is responsible for the absorption of the laser radiation in plasma. 
This condition is failed, if the influence of the strong EM wave lead to development of an 
instability. Hence, pulse duration r of an EM wave has to satisfy the condition fir ^ 1, 
where /i is the maximal increment of the instability of the plasma in the strong laser 
held. 

In general, analytical integration over momentum po and scattering phase ips is 
impossible, and one should make numerical integration. The latter along with derivation 
of asymptotic formulas for absorption coefficient a will be done in the next section. 

3. Numerical Treatment: Asymptotic Formulas 

As we are interested in superintense laser pulses of relativistic intensities, then it is 



convenient to represent the absorption coefficient (26) in the form of dimensionless 
quantities: 

= o^n ^^2\f2 / ^Po/(Po) / — -=3 

"0 27r(l + C^)^^ J Jo p3(^,)A 
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Here 

ao = 4ZVf A^n^ne, (32) 



and Te is the classical electron radius. In Eq. (31 ) the dimensionless momentum, energy, 
and temperature were introduced as follows: 



e 

mc mc^ mc^ 



and the dimensionless relativistic intensity parameters of EM wave, 

^ii^s) = ^o(eicos?/'s + e2Csin?/'^). 
The scaled relativistic distribution function is 

7o 



and Coulomb logarithm 

where Lq = In (A/27rre2')^. In the near-infrared and optical domain of frequencies, at 
the Z=l — 10, Lo~30. Taking into account that normalized absorption coefficient 



(31) depends on Z and u through logarithm Lq, for the numerical simulations we will 
not concretize Z and u, assuming Lq ~ 30. 

It is well known that the kinematics of an electron in the field of a strong EM 



wave essentially depends on the polarization of the wave 10 . In particular, for the 
particle initially at rest, in the CPW, energy 7(V's) and momentum p [ips) are constants, 
since ^^(^/'s) = const. Meanwhile in the LPW p{ips) oscillates and as a consequence 



small values of p{'ips) give the main contribution in Eq. (31). The latter leads to more 
complicated behavior of the dynamics of SB at the linear polarization of a stimulating 
strong wave. Besides in the case of CPW, thanks to azimuthal symmetry one can make 
a step forward in analytical calculation and obtain explicit formula for the absorption 
coefficient at superstrong laser fields. Thus, taking into account azimuthal symmetry in 
the case of CPW, one can make integration over phase ips, which results in 



ao 2 



j dPof (Po) (^7 (Po) ^ + 1^ ^0 



where 



X ^^Le., (34) 



7 (Po) = 7o + i [^0 + 2Po^o sin cos <^o] , (35) 
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Figure 2. (Color online) Total scaled rate of inverse-bremsstrahlung absorption (in 
arbitrary units) of circularly polarized laser radiation in Maxwellian plasma versus the 
dimensionless relativistic invariant parameter of the wave intensity for various plasma 
temperatures. The wave is assumed to be circularly polarized. 



At the large >> 1, taking into account Eqs. (35) and (36), from Eq. (34) one can 
obtain 



a 



so 



dVof (Po) 



A 



cb ' 



(37) 



where 



-(c) 
^cb 



so 
2A 



The formula (37) shows the suppression of the SB rate with increase of the wave 
intensity. Ignoring weak logarithmic dependence, we see that absorption coefficient 
inversely proportional to laser intensity: a ~ l/^o- -^^^ large the dependence of 
the absorption coefficient on temperature comes from A in Eq. (37). In particular, for 
initially nonrelativistic plasma Tn << 1 in Eq. (37) one can put A ~ 1, which gives 

2 



a = ac 



so 



In 



(3J 



The relation for the absorption coefficient in the case of LPW is complicated and even 
for large one can not integrate it analytically. Therefore, for the analysis we have 
performed numerical investigations, making also analytic interpolation. 

The results of numerical investigations of Eq. (31) are illustrated in Figures 2-7 
both for the CPW and LPW. 

To show the dependence of the inverse-bremsstrahlung absorption rate on the laser 
radiation intensity in Fig. 2 it is shown scaled rate a/ao versus relativistic invariant 
parameter of the wave intensity for various plasma temperatures. The wave is assumed 
to be circularly polarized. As is seen from this figure the SB rate is suppressed with 
increase of the wave intensity and for the large values of .^o it exhibits a tenuous 
dependence on the plasma temperature. The behavior is also seen from Fig. 3, where 
total scaled rate of inverse-bremsstrahlung absorption of CPW in plasma, as a function 
of the plasma temperature is shown for various wave intensities. Here for large 
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Figure 3. (Color online) Total scaled rate of inverse-bremsstrahlung absorption (in 
arbitrary units) of circularly polarized laser radiation in plasma, as a function of the 
plasma temperature (in units of an electron rest energy mc^) is shown for various wave 
intensities. 




Figure 4. (Color online) Density plot of the total rate of inverse-bremsstrahlung 
absorption scaled to asymptotic rate ac (in arbitrary units), as a function of the 
plasma temperature (in units of an electron rest energy mc^) and the dimensionless 
relativistic invariant parameter of the circularly polarized laser beam. 




Figure 5. (Color online) Total scaled rate of inverse-bremsstrahlung absorption (in 
arbitrary units) of linearly polarized laser radiation in plasma versus the dimensionless 
relativistic invariant parameter of wave intensity for various plasma temperatures. 
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Figure 6. (Color online) Total scaled rate of inverse-bremsstrahlung absorption (in 
arbitrary units), as a function of the plasma temperature (in units of an electron rest 
energy mc^) is shown for various wave intensities. The wave is assumed to be linearly 
polarized. 
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Figure 7. (Color online) Density plot of the total rate of inverse-bremsstrahlung 
absorption scaled to asymptotic rate ul (in arbitrary units), as a function of the 
plasma temperature (in units of an electron rest energy mc^) and the dimcnsionless 
relativistic invariant parameter of the linearly polarized laser beam. 



values of we have a weak dependence on temperature, which is a result of laser 
modified relativistic scattering of electrons irrespective of the initial state of electrons. 
The absorption coefficient a decreases as 1/^q in accordance with analytical result (37). 
To clarify the range of applicability of the asymptotic formula (38) in Fig. 4 density 
plot of the total rate of inverse-bremsstrahlung absorption scaled to asymptotic rate ac, 
as a function of the plasma temperature and the relativistic invariant parameter is 
shown for CPW. As is seen in the wide range of Tn and one can apply asymptotic 
formula (38). 

In Fig. 5 and 6 it is shown the total scaled rate of inverse-bremsstrahlung absorption 
of LPW in plasma versus the dimensionless relativistic invariant parameter of wave 
intensity and temperature, respectively. With the interpolation we have seen that a 
decreases as l/^o''^ exhibits a tenuous dependence on the plasma temperature. 
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Making analogy with the case of CPW, for the large we interpolate a by the following 
formula: 



a 



ftp 

9f5/4 



In 



(39) 



As is seen from Fig. 7, in the case of LPW and for the moderate temperatures, with 
the well enough accuracy one can apply asymptotic rate (39). 



4. Conclusion 



We have presented a theory of inverse-bremsstrahlung absorption of an intense 
laser radiation in relativistic Maxwellian plasma in the relativistic low-frequency 
approximation. The coefficient of nonlinear inverse-bremsstrahlung absorption has been 
calculated for relativistic Maxwellian plasma. The simple analytical formulae have been 
obtained for absorption coefficient at asymptotically large values of laser fields both for 
circularly and linearly polarized radiations. The obtained results demonstrate that the 
SB rate is suppressed with the increase of the wave intensity and for large values of 
absorption coefficient a decreases as for circularly and as for the linearly 

polarized one in contrast to nonrelativistic case where one has a dependence 1/^q (Tt) . 
The SB rate is suppressed with increase of the plasma temperature but for the relativistic 
laser intensities it exhibits a tenuous dependence on the plasma temperature. 
This work was supported by SCS of RA. 
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